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ABSTRACT 

One might expect light to be scattered when it passes through a gravitational wave, and might 
hope that in favourable circumstances these scatterings could be observed on Earth even if 
the interaction occurs far away. Damour and Esposito-Farese, and Kopeikin, Schafer, Gwinn 
and Eubanks, found that there were cancellations making such effects disappointingly small. 
Here I show that those cancellations depend on the emission of the light occurring far behind 
the gravity-wave source; for light-emissions near that source, larger effects are possible. I first 
develop a covariant treatment of the problem in exact general relativity (the propagation of 
light being modelled by geometric optics), and then specialise to linearised gravity. The most 
promising candidates identified here for detection in the not-too-distant future would involve 
sufficiently tight binaries as sources of gravitational radiation, and nearby pulsars as light- 
sources. In some favourable but not extreme cases, I find offsets in the pulses' times of arrival 
at Earth by ^ 10^^" — 10^^ s, with periods half the binaries' periods. 
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1 INTRODUCTION 

The detection of gravitational waves is a critical test of Einstein's 
theory of general relativity. We have at present good indirect evi- 
dence for the existence of such waves in the orbital decays of binary 
pulsars (Taylor & Weisberg|1989^ , but no direct evidence. Unfor- 
tunately, we expect most astrophysical sources to be so distant that 
the waves should be very weak by the time they reach the solar 
system, making terrestrial detection a challenge. 

It is natural to ask whether gravitational waves might have ef- 
fects on matter or light closer to the sources, where the waves are 
stronger, but with the effects being observable on the Earth. In par- 
ticular, one might search for modulations in light caused by its pas- 
sage through gravitational waves, for instance light from another 
electromagnetic source passing nearby a neutron star or binary sys- 
tem. (Here and throughout, 'light' means electromagnetic signals 
of any frequency which can be treated by geometric optics.) These 
modulations could include variations in intensity, in apparent po- 
sition, in polarisation, in phase, or in time-dilation (often referred 
to in the literature as time-delay). Of course, such variations could 
also be produced by other means; one would need further informa- 
tion to conclude they were caused by gravitational waves. 



1.1 Previous work and its implications 

Such possibilities were contemplated by a number of authors, no- 
tably |Sazhin| ( |1978) , |Labeyrie| p993) , |Durrer| ([1994) and [F^ 
( |1994 1; it seemed from these papers that some of the effects might 



be detectable in the near future. Concerns, though, were raised 



It^ 

'99]l 



E-mail: helfera@missouri.edu 



about the adequacy of those treatments, and Damour & Espositi 
Farese (1998) and Kopeikin, Schafer, Gwmn & Eubanks (1999] 
(see also ^Kaiser & Iaffe^^l997^) , after more considered analyses, 
found very much smaller effects than then the optimistic ones in 
the earlier papers]^ The effects were so small that they seemed to 
preclude the detection of gravitational waves by these means, at 
least for known classes of systems in the near future. This conclu- 
sion has been accepted by much of the community ( |Schutz|2010[ 
[Book & Flanagan|2011| ). 

But the situation is not so clear-cut. For one thing, the later 
authors' investigations, while overlapping the earlier ones, do not 
cover all of the physical possibilities they considered. Also the re- 
sults as they stand strongly suggest there are essential points which 
have not yet been understood. 

We can appreciate this if we ask how such gross discrepancies 
between the earlier and later work could have occurred. While it is 
true that much of the earlier work involved rough approximations, 
one would think that in order to explain such a large and systematic 
change there must have been a shift in what was taken to be the 
essential physics. 

The earlier papers typically and very plausibly found that 
the light rays responded to the gravitational radiation fields of 
the sources, but later authors found that there were cancellations 



^ The literature on this subject is extensive; for more references and dis- 
cussion, see IKaiser & Jaffel (TW/J; IDamour & Esposito-Faresel JT998); 
|Kopeikin et al. 1 1999i. Additional references, on the light's polarisation, 
are, Kopeikin & Mashhoo n, ^2002j ; ,Prasanna & Mohanty,^2002j ; |Faraoni] 
|2008|; andlLesovik et al.|j2005| on the Ught's phase. 



© 2012 RAS 



2 Adam D. Heifer 



which resulted in those terms not entering [J Indeed, Damour and 
Esposito-Farese pointed out that pure, vacuum, gravitational waves 
would give no overall scattering in linearized theory. More specif- 
ically, both Damour and Esposito-Farese, and Kopeikin et al., re- 
ported that what scattering of a light-ray by an isolated source (on 
a straight world-line in linearised theory) does occur depends only 
on the source's configuration at the light-ray's point of closest ap- 
proach. This is something which cannot be locally causally deter- 
mined, and therefore must depend on a non-local cancellation of 
effects over the light-ray's trajectory. 

These observations suggest that the cancellations are tied to 
the infinite extent of the light-ray's trajectory, and may not occur for 
propagation along finite, or only half-infinite, segments|^We shall 
see below that this is indeed the case. In other words, the cancella- 
tions are present when the light effectively comes in from infinity, 
but larger effects are possible when the emission is from a finite 
point, in particular when it might be close to the gravitational-wave 
source. Such situations are very much of interest, because there 
are many astrophysical configurations in which there is a source of 
light near an expected source of gravitational radiation. 

This fits well with other points in the literature. One can see 
that many of the earlier works did, either explicitly or in effect, 
use scattering estimates which depended on integration over half- 
infinite or finite segments rather than doubly-infinite lines, and in 
some cases (depending on the situation considered) the authors 
were correct in wanting to do so. And Damour and Esposito-Farese 
(see also |Faraoni|2008 > correctly noted that overcoming the cancel- 
lations would depend on 'edge effects'. 

The paper of Kopeikin et al. is sometimes considered to show 
that no significantly larger effects arise from finite or half-infinite 
segments. This, however, is not really correct. While these authors 
did indeed work out formulas for the null geodesies around cer- 
tain sources in linearised gravity, those formulas are very lengthy 
and were examined in detail only in certain limiting cases. None 
of those cases corresponds to the configuration of ultimate interest 
here, a light-emission reasonably close to the gravitational-wave 
source (precisely, within of the order of the impact parameter to the 
light-ray's point of closest approach to the gravity-wave source) re- 
ceived by a distant observer. In fact, the computations of Kopeikin 
et al. are consistent with the possibility of larger effects from re- 
stricted segments. 

There is, additionally, a key issue which must be properly ad- 
dressed if we are to take into account the emission (or reception) 
of light-rays at finite events, especially events near gravitationally 
radiating sources. We do not expect to learn anything about grav- 
itational radiation from the reception of one light-ray; rather, it is 
in the changes in the light received over time that the we hope the 
information is encoded. In other words, it is differential scattering 
effects, of successive light-rays relative to each other, which are of 
interest. But to treat these properly, we must know we are mod- 
elling the emissions of the successive rays accurately. This is po- 
tentially problematic when the light-sources themselves are in the 
region affected by the gravitational waves, because the effects of 

^ It is worth noting that essentially the same feature comes up in the text- 
book analysis of the geodesic deviation of linearised plane waves in a stan- 
dard gauge, used to discuss gravitational wave detectors, where one finds 
that the deviation depends on the metric coefficients at the end-points only 
of the geodesic (Misner et al. 1973 1. 

^ Since 'light-ray' has an established meaning as a path, whether finite, 
half-infinite or doubly-infinite, taken by light in geometric optics, I will not 
use 'ray' to mean a half-infinite segment. 



the waves on the light through their perturbations of the motion of 
the light-sources compete with the other contributions. 

Another way of saying this is that the received signal depends, 
not only on the gravitational field through which the light has 
passed, but also on the apparent acceleration of the light-emitter, as 
inferred from electromagnetic measurements at the receiver. Even 
if the emitter suffers no local acceleration, its trajectory as recon- 
structed from light-signals coming into the detector will appear 
to involve accelerations, because of the gravitational radiation the 
light has passed through. (And so one might want to say that the 
variations in the received signal are not purely due to scattering, 
but also to this non-local acceleration of the source. However, from 
a general-relativistic point of view these are two facets of the same 
thing.) 

This is a point which has only been partially attended to in the 
literature on scatteringj^In certain restricted cases, where the grav- 
itational waves have a simple form and the light-rays considered 
are aligned to take advantage of it, the issue has been considered 
(e.g. ! Sazhin|[l978[ l, but the pattern has been to compute the null 
geodesies but not examine all the details of the time-dependence of 
light-emission and reception. 

For instance, [Kopeikin et al.| l |1999[ l worked out the null 
geodesies around certain gravitational sources in linearised the- 
ory, in terms of the geodesies' initial positions and tangents. They 
derived formulas for the coordinates of finite points along the 
geodesies, but they did not fully resolve the question of how to 
choose the data for a one-parameter family of null geodesies emit- 
ted at one world-line and received at another. (They addressed 
some aspects of the issue, with their discussion of the appropriate- 
ness of their gauge choice. That argument shows that some light- 
emitters may be modelled by taking the spatial coordinates, in the 
gauge used, fixed. For such an emitter, the initial positions of the 
null geodesies would simply be the points on the emitter's world- 
line. However, this leaves open the question of how to choose the 
geodesies' tangents, so as to ensure the geodesies go from the emit- 
ter to the receiver; it also leaves open the question of how to treat 
generic emitters, which need not have their space coordinates con- 
stant.) Developments of this approach can of course take into such 
variations; see Kopei kin et al.| ( (20Il[ l. 

We wish then to extend the existing treatment of light- 
scattering by gravitational fields to propagation over finite or half- 
infinite segments, accounting for the gravitational waves' effects on 
the light-sources in those cases. It is worthwhile considering what 
the best approach is. 

1.2 A more general setting 

Virtually all work that has been done in this area has been done 
in the linearised-gravity (or at least linearised gravitational-wave) 
setting, and has used detailed coordinate calculations in specific 
gauges. There are reasons for wanting a more general, and invari- 
ant, approach. 

Most importantly, the core issue - how light emitted from suc- 
cessive events along one world-line is received along another, after 

Curiously, there is another line of thought on gravitational-wave de- 
tection, where the proper resolution of this point is an integral, if little- 
emphasized, part of the foundations. This is the idea of pulsar timing ar- 
rays. Often the basic principle is explained as the use of the pulsar and the 
Earth as the two ends of a gravitational-wave detector: thus the emissions of 
pulses ai'e from events along the pulsar's world-line, which is freely falling 
(e.g. |HeUings & Downs] l"983} . 



© 2012 RAS, MNRAS OOO.pTfIS] 



Light, gravity and pulse-time offsets 3 



passage through a region of space-time - is evidently a basic prob- 
lem in relativity: it is worth treating in generality. It is a tenet of 
relativity that the physics is invariant, so invariant treatments focus 
most directly on the physics; a central goal of this paper is to give 
one. 

As emphasized earlier, the key physical observables are asso- 
ciated with differential scattering, of successive light-rays relative 
to each other. This is governed by the geodesic deviation equation, 
a linear ordinary differential equation involving the curvature, and 
it is the solution of this equation which is the central mathemat- 
ical chore in an invariant approach. By contrast, gauge-based ap- 
proaches focus on the metric and on integrating the null geodesies, 
a considerable task. It is only once these are known that one can 
proceed to extract the measurable differential effects (and, as noted 
above, this further step has not generally been completely worked 
out in these approaches). One sees that the invariant approach goes 
more directly to the observables. Also, by limiting oneself to dif- 
ferential effects, one expands considerably the cases which can be 
treated, since one need not confront the problem of finding the null 
geodesies explicitly. The invariant formulas derived here would, for 
instance, be applicable in a cosmological context, or to treat light- 
rays in strong gravitational fields. 

With the invariant formulas in hand, one can always specialise 
to gauge choices which may be convenient for particular prob- 
lems. But not having to make such a choice too early saves much 
work, and does away with many conceptual issues associated with 
gauges. 

In some of the literature the motivation for the choice of gauge 
is to be confident that one is indeed modelling gravitational waves, 
as opposed to other gravitational disturbances. In this paper, that is 
a distinction which is best made, not at the beginning of the analy- 
sis, but rather once we have the general formulas. This is because 
there is no very simple invariant characterization of what counts as 
gravitational radiation and what does not. We first derive general, 
invariant, formulas; we may then specialise these to systems which 
are gravitationally radiating. 

1.3 Goals and results 

The aims and results of this paper are: 

(a) To provide an invariant treatment of the effects of gravity 
on light, in the approximation where the electromagnetic field can 
be treated by geometric optics, in exact general relativity, and a 
specialisation of that treatment to linearised gravity. In particular, 
while much of the motivation for this work comes from wanting to 
understand gravitational waves, the approach here allows us to treat 
any gravitational fields in general relativity. 

I emphasized above that the effects of interest were all dif- 
ferential ones, where one light-ray is compared with a nearby one. 
Thus the main tool in the analysis is the geodesic deviation equa- 
tion. Essentially, the task is to solve this equation with initial and 
final conditions corresponding to the emission of light at one world- 
line and its absorption at another. 

This approach allows us to disregard many of the coordinate 
issues which occur elsewhere. For instance, previous papers have 
tended to discuss 'time-delay' effects in terms of local coordinate 
times, whereas it is really the clock times of the emitter and the 
receiver of light which are relevant; the approach here gives an 
integro-differential equation directly relating the clock times. 

There are similar equations for the changes observed at the 
receiver in the wave-vector (this is equivalent to equations for the 
change in phase, since the wave-vector is the gradient of the phase) 



and the electromagnetic field (and thus the amplitudes and polari- 
sations). In all of these cases, the changes are determined partly by 
the geometry of the space-time the light propagates through, partly 
by the world-lines of the emitter and receiver, and partly by any in- 
trinsic changes in the field at the emitter. All of these contributions 
are accounted for directly. 

These results are closely related to those of Damour and 
Esposito-Farese; the chief difference is that they worked with the 
field's Fourier transform and gave global results, whereas our anal- 
ysis is directly in terms of the radiative data and considers not only 
the limiting case of a doubly-infinite geodesic, but how that limit 
arises. Similarly, the results here are also compatible with the com- 
putations of Kopeikin et al. 

(b) To give rough estimates of the effects in question in the 
case of quadrupolar radiation in linearised gravity. 

Of the various effects on the propagation of light by its passage 
through, or emission in, a gravitational wave field - changes in its 
position on the celestial sphere of the receiver, in its phase, polari- 
sation, amplitude, or time-dilation - the last is, at present, the least 
difficult to imagine measuring. We shall find that changes in the 
position on the celestial sphere and phase are suppressed by a geo- 
metric factor (related to the 'lighthouse effect'). While changes in 
polarisation or (fractional) amplitude might be, in numerical terms, 
of the same magnitudes as those due to time-dilation, it is hard to 
foresee polarisation or intensity measurements accurate enough to 
reveal gravitational- wave effects. 

The time-dilation effects (often referred to as 'time-delay ef- 
fects') arise from the distortion ri(r) of the light-emitter's proper 
time Ti as measured by light-signals arriving at the receiver's 
proper time r. This function ri (r) depends on gravitational radi- 
ation the light encounters; in particular, we shall see that an in- 
tegral of a curvature component over the geodesic contributes to 

One result of this is a red-shift; also, if the emitter can be ex- 
pected to give off signals at known proper times (for instance, if it 
is a pulsar), then the light's passage through a gravitational wave 
may cause the times of arrival of those signals to wander slightly 
from what they would have been, had no radiation been present. 

The detection of these effects would require training a tele- 
scope, whether optical or radio, on a light-emitter in the vicinity 
of a gravitational-wave source. For this reason we are interested 
primarily in sources which have at least a reasonable likelihood of 
radiating in a reasonable observing period. These will be continu- 
ous, localised, sources of gravitational waves, which unfortunately 
are thought to be weak. 

The best sort of candidate found here would involve a tight 
binary as the source of gravitational waves, with a nearby pulsar 
(perhaps as a distant tertiary) as the light-emitter. For two 5Mq 
stars (or black holes) in a .1-day circular orbit, I find possible dis- 
placements Armax in pulse arrival times ~ 2 x 10~^° s. For (as 
might possibly be found in a globular cluster) a one-solar-mass star 
in a one-day circular orbit around a 10'' Af© black hole, one has 
Armax ~ 3 X 10~' s. For (as might happen near the galactic centre) 
a star of mass m in a ten-year circular orbit around a 4 x 10^ Mq 
black-hole, one has Armax ~ 7 x lO~^°(m/M0) s. These dis- 
placements would vary with a period half the orbital period. 

1.4 Organization 

The next section develops a general formalism for the analysis, 
valid in full general relativity, based on geodesic deviation. In Sec- 
tion 3, this is specialised to the case of linearised gravity; there the 
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geodesic deviation equation can be solved in terms of curvature 
integrals. Section 4 considers the effects of quadrupole sources (al- 
lowing both 'electric' and 'magnetic' terms) of gravitational radia- 
tion on light-rays in the sources' wave zones. In Section 5, those re- 
sults are applied to estimate the magnitudes of the effects for some 
expected sources of gravitational radiation. Section 6 contains a 
summary and discussion. 

The reader wishing to skip the derivations will find the for- 
mula (jTTJ for the comparison of emitter and receiver clock-times at 
the end of section 4; the meanings of the symbols in this formula 
are recapitulated in the paragraph containing it. The rough estimate 
of the magnitude of the effect is \79\ . Sections 5 and 6, dealing 
with the estimates, depend (except for a few technical comments 
which can be skipped) only on this rough form and can be read 
independently of the earlier sections. 

Notation and conventions. Except where otherwise specified, 
the notation and conventions here are those of [Penrose & Rindlerl 
( |1984[[T986[ l, which will also serve as a reference for all material 

not otherwise explained. The metric signature is H , and the 

curvature satisfies [Va, Vb]v'' — Rabc^'v'^ ■ Factors of c, the speed 
of light, are omitted until the end of Section 4. 

In keeping with this paper's main aim of an invariant treat- 
ment, the basic formalism is, in principle, the abstract-index one 
of Penrose and Rindlen However, in fact, in almost all cases the 
tensorial indices in this paper may be interpreted equally well as 
abstract indices or as component indices with respect to a chart. 
Where there is any difference, it is noted. 



2 GENERAL FORMALISM 

Let us suppose we have two world-lines '^jijj) (with j = 1, 2) in 
space-time, with Tj proper time on each. These will be the world- 
lines of the source and the detector of the light-rays [^They need not 
be geodesic, and they need not be in any asymptotic or weak-field 
region. Let p(s, r) be a smooth family of light rays from 71 to 72: 
for each fixed r, the light ray runs from 7i(Ti(r)) = p{si{t), t) 
to 72(T"2(r)) — p{s2{t), t), with s an affine parameter along the 
ray. Then 1°" = dsP will be the tangent null vector, and it)" — drP 
will be the connecting, Jacobi, fieldj^Note that we have 



d7 



7i 



72 



dS2 

d7 



r + 1 



(1) 

(2) 



on 71, 72. We shall take r = r2, so as to index the light-rays by 
the receiver's proper time. There is some normalisation freedom in 
this: for each r, the vector (s2(r) — s\{t))1"' is determined, but 
the individual values S2 (r), si (r), are not. It will be simplest to 
take 



Si(r) = const , S2(t) = const , 



so 

dri 
dr2 



7x = lu on 71 



and 



72 = If on 72 . 



(3) 



(4) 



^ We should write 7° (tj ) for the coordinates of a curve in a chart, with 
a a coordinate index; similarly we would have p"(s, r) for the coordinates 
of a family of light-rays below. (Since the coordinates themselves are not 
tensorial quantities, the quantities 7j (tj ) and p(s, t) do not carry abstract 
indices. 

^ So these relations would be = flsp" and w'^ = drP"^ in a coordinate 
chart. 



We shall be interested in how the light-ray varies with r, 
and hence in the connecting field w°'. Contracting the equation 
/ • Vw" = w • Vr with we find / • V(it) -0 = 0, which, 
with eq. |4](, gives us 



dri 
d7 



/ ■ 71 = / ■ 72 . 



(5) 



It is dri / dr which gives rise to what are often called 'time-delay' 
effects, and eq. l[5| will allow us to solve for these. 

We shall want the connecting vector w"" in terms of the curva- 
ture along the light-rays. For this, we must use the Jacobi equa- 
tion in some detail. Let V^tiSyT), yf,(s,r) be solutions, so 
(l ■ \/fU\ = Fl''Rpcg''U''b, {I ■ V) V% = Fimpcg^V^b, with 



C/''b(si,r) = (5"i, 
/•VI7%(si,r) = 0, 

and 

V\{s^,r) = 
/•VV\(si,r) = 

Then the connecting field is 

w =U bCt +V bP , 



0, 

5\. 



(6) 
(7) 



(8) 
(9) 



(10) 



for some a", /3" (elements of the tangent space at 71 (ri)). We can 
find Of", Z?" by using eq. l|4|l; this gives 



pa 



{V{S2) 



■ b 

72 - 



-^U{S2) e7l 



(11) 

(12) 



Note that I ■ V{1 ■ w) = 0, applied to eq. 1 lOi and evaluated at 
s = s-i, implies I ■ jS = Q. 

Of course, at a conjugate point V^b will not be invertible; con- 
jugate points will be discussed elsewhere. 



2.1 Time dilation 

Much of the literature is phrased in terms of 'time delays' , the delay 
being taken to be the difference in the coordinate times of emission 
and reception. This is not an invariant concept, and so it is then 
corrected (or corrections are at least as a matter of principle con- 
sidered) to take into account differences between the clock times 
and the coordinate timesQ 

We will work with the clock times directly, the basic quantity 
of physical interest being ri (r), the time along 71 at which a signal 
was sent to arrive at time r on 72. Then dri/dr, the time-dilation 
or red-shift factor[^It is changes in this quantity which may carry 
the imprint of gravitational radiation. 



^ The literature also refers to these as 'Shapiro time delays'. But in the 
Shapiro case the configuration is different: the radio signal is reflected off 
a planet and then received by a detector which has (essentially) the same 
world-line as the emitter, and so there is a well-defined time-delay between 
emission and reception. It is this difference in clock times which Sha piro| 
(19641 originally called the time delay. The present techniques can easily 
be adapted to give an invariant formula for that. 

* One has dri /dr = {1 + z)~^, where z is the usual measure of red-shift 
of signals from 71 to 72 . 
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We have 



I ■ 71 

= (1 ■ 7i)~^ {1 ■ 7iw • v(; • 72) - r 72m ■ v(; ■ 71)) 

= (l ■ 7l)~^ (' • 7l72a^ ■ Vw" - I ■ 727lai ■ Vto") 

+ {l ■ 7i)~^^ ■ 72 - ■ 7i)~^(^ • 'yifl ■ 71 
= {I ■ i {wj -Vw")} — {wj -Vw")] 

\ Is — S2 

+ {l ■ 7i)~^^ ■ 72 - (/ • 7i)~^(^ • 72)^^ ■ 71 
Expressing this in terms of U and V, we have 



(13) 



(/ ■ 7i)''72ar V((7%a'' + V%/3') 

I s — S2 

-(/ ■ 7i)~^(/ ■ 72)71 • P 

+ (1 ■ 7i)~^^ ■ 72 - 7i)~^(^ • 72)^^ ■ 71 ■ 



(14) 



The last terms of course vanish when the source and the emitter are 
freely falling. (In eqs. l |13[ ). and subsequently, the accelera- 
tions 'jj ~ "fj ■ V7" are taken with respect to the proper time Tj 
along the corresponding world-line.) 



2.2 Change in wave-vector 

In the geometric-optics approximation, the wave-vector is 

fc''=a;ir/r7i, (15) 

where uji is the angular frequency with respect to the frame of the 
emitter. We are interested in the time-dependence of fc" at the de- 
tector's world-line, 72(1"). For simplicity, we will assume here that 
uji is constant. (Otherwise, in what follows, one simply gets an ex- 
tra term, from the product rule.) Then applying to ■ V to eq. jl5[ l, 
we have 

72 • Vfc" = Ldl{l ■ 7l)~^ {I ■ 'jlW ■ Vi" — l"'W ■ V(i ■ 71)) 

= ■ 7i)~^ I I ■ ■yil ■ Vto" 

—l"^ {w ■ V lb)^\\ — r/fcif • V71 



= a;i(Z ■ 71)"' {I ■ ■ V{U\a' + V\p')\ 

-r [/?.7i + (r^i)-i(r^2);,7j] ) . (i6) 

The last term in the brackets is proportional to the acceleration of 
71 , and vanishes if that world-line is freely falling. Also, the equa- 
tion l|14[ ( for the time-dilation can be regarded as a consequence of 
eq. l[f6K since d^ri/dr"^ = 72 ■ V(72 • k/uji). 

Finally, a physicist receiving signals must decide how to com- 
pare successive measurements of fc" along 72. If 72 is a geodesic, 
there is a natural choice: parallel- transport, which leads to the dif- 
ferential formula (jT6|. If, however, 72 is not a geodesic, one might 
prefer to use Fermi-Walker transport. It is easy enough to inter- 
convert the two, the Fermi-Walker derivative of k°' along 72 being 



derivative gives us ones less sensitive to the geometry of the world- 
line 72 ; the Fermi- Walker derivative is more frame-dependent but 
measures more directly the changes in frequency and spatial wave- 
vector relative to the observer. The same principle will apply to 
changes in the field. 



2.3 Change in the field 

We may also analyse the effects on the received electromagnetic 
field of its propagation through the gravitational field; this includes 
changes in amplitude and polarisation. 

As is well-known, in geometric optics the field Fab is trans- 
verse to the direction of propagation. That is, however, not a rel- 
ativistic formulation of the condition; relativistically the statement 
is that must be a repeated principle null direction of Fab ( |Pirani| 
[T965] [Penrose & Rindler|1986) . We may express this conveniently 
by choosing a complex null vector m" (covariantly constant along 
and normalised to i ■ m = 0, m ■ m = —1); then the field must 
have the form 



Fab = (t>{lamb — rrialb) + conjugate 



(18) 



for a scalar field (j). One can fix the freedom in the choice of m" to 
have </> ^ 0; then the two-form ZaWi6 — ?7ia/f, carries the polarisation 
information. One has 



(19) 



where p = — (1/2)V ■ I is the convergence of 1°'. 

It is convenient to introduce a luminosity distance v such that 
l-Vx = — px. Then I ■ V {xFab ) ~ and xFab is parallel-transported 
along the null geodesic. This means that the field Fab diverges as 
1/x at 71, but that is simply a mathematical artefact of modelling 
the emitter as a point source. Really, one should imagine a finite 
surface of emission, and xFab taking direction-dependent limits as 
one approaches this surface. 

The natural normalisation for r is then with respect to the 
world-line of the electromagnetic source, that is (7i-/)~^/-Vr — 

1. The direction-dependent limit of xFab on 71 is a measure of the 
intrinsic strength of the field at the source, which will vary along 
the world-line. 

We can express r in terms of the quantities already given. 
Since p measures —1/2 the logarithmic rate of increase of the sur- 
face area element along the rays abreast the luminosity distance 
is the square root of the area element. Taking into account the nor- 
malisation, we have 



■ = (71 • l)V -'iV"- cV'' dmafribml'^fn'^^ 



(20) 



72 • Vfc" + (7272 - 7272)^6 , (17) 

so one would supplement l |16| ) by an additional term. The choice 
of which quantity to use is really a question of which measures 
of the change one is most interested in reporting. The covariant 



We have then 
l-VwV (xFab) 



= [l-\/,W-\/]xFab 
= -2Fw''Rpq[a'''<:Fb]c 



(21) 
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Integrating this, we have 

72 • V (vFab) = W-V (vFab) 
= W-V(xFab)\ 



fS2 

_2 / lP-,n1 



l''w''Rpq[a''xFt]cds 



d7 



71 • V {vFa. 



-2 / Fw''Rpg[a''^Ff,]^ds 



d7 



71 ■ V {vFab 



Pw^Rp^u ds tFt 



(22) 



where we understand that parallel propagation along the null 
geodesic has been used to identify w ■ V{tFab)\ with a tensor 

Eisi 
Then 



72 • VFab 



. _ , jl-\/{xFab)\ 

72 ■ Vr , dn Ui 

^ab H ; 

r dr r 



Fw''Rpgla'' ds ) Fftjc . 



(23) 



The last term is a pleasingly clean Lorentz transformation derived 
from a curvature integral; the middle term corresponds to changes 
in the received field due to intrinsic changes in the source. The 
first term, due to changes in luminosity distance, can be expressed 
(somewhat lengthily) in terms of the data we have, as follows. 
We have, from eq. l |20| >, 

72 • V log r 



Here 

72 • V71 • I 



72 • V log(7i ■ (24) 
+ (1/2)72 ■ V log V°cV''dmam6m'''m'*' . 

w ■ V71 ■ I 

(w ■ \/ji)la + jlW ■ Via 



dri „ 



dr 



7i • ; + 7i • /3 . 



(25) 



Differentiating the normalisation conditions for m", we find 

72 • Vm^Jfib] = w ■ Vm^Jfib] 

= -(/• VC/^'^a' + Z- VV%/3«)mpn[,m6j 
— conjugate 

+ terms proportional to la , h , (26) 

where the terms proportional to la or lb will not contribute to eq. 
l |24| , because those are eigenvectors of Vq. Finally, we must com- 
pute 72 • Vl^";, — w-W^bl • Differentiating the Jacobi equation 
(which V^b satisfies) and working out some commutators, we find 

I ■ V^w • VV"b - Fl^Rpcq^-w ■ VV'b = 5"(, (27) 



^ Thus an integral written as f"^Qab{s)ds is really 
f^^ Pa''{s)Pb'''{s)Qcd{s) ds, where Pa'^{s)\c is the result of parallel- 
propagating Ac to 72 along the null geodesic p(s, r); similarly we should 
haveP„<=(si)P6<*(si)M,- V(rF,d)| for «, ■ V(rF„i,) I . 



where 

S\ = {w-V{Fl''Rj,,q''))V''b-w''l''Rpq,''l-VV\ 

-I ■ Viw^l^Rpqc^Vb) . (28) 

We can regard eq. \H) as a sort of inhomogeneous Jacobi equation 
with source S°'b, and solve it by variation of parameters. The result 
is 

wVV = [ [U{s)~V{s){y'^){s)U{s)\x (29) 

J so 

[l-VU ~{l- W)V-^U] (s)S'(s) ds , 

where the indices have been omitted (and matrix operations are to 
be understood throughout) in the interest of clarity. 

Combining eqs. \25) , \26) , \29) (evaluated at s = S2) and l |28[ > 
gives the first term on the right in eq. l |23| l, and this completes the 
formula for the change in the field at the receiver. 



3 PASSAGE TO LINEARISED GRAVITY 

We will now specialise to linearised gravity. Thus we regard the 
metric as a first-order perturbation of the Minkowskian one. We 
shall nevertheless avoid an explicit choice of gauge, and continue 
to present the results in an invariant form, in order to keep the ge- 
ometry and physics as clear as possible. 

There is an overall issue to keep in mind in such schemes: 
in general, a quantity of interest will have both zeroth- and first- 
order terms. If it is not a scalar, and its zeroth-order term is non- 
vanishing, then a first-order gauge change will in general add in a 
portion of the zeroth-order term to the first-order term. This means 
that the decomposition into zeroth- and first-order terms is not in- 
variant unless the zeroth-order term vanishes. 

In what follows, many of the contributions we compute will 
be purely first-order, and thus will have invariant interpretations. 
However, for each of the effects there will also be zeroth-order 
terms. (For instance, if 71 and 72 are skew time-like geodesies in 
Minkowski space, then ri(r) will incorporate a time-dependent 
Doppler effect.) Thus there will be certain contributions which 
have no invariant decomposition into zeroth- and first-order effects; 
for these, any attempt to specify them in terms of a background 
Minkowski geometry will require choosing a gauge. 

The place this mixing of zeroth- and first-order terms will 
show up is when we use parallel transport along the null geodesic 
to identify the tangent 7" to 71 with a vector at 72. In fact, the 
quantity which will enter is 



5" 



72 



dri a ■ b 
d7^ "^^ ■ 



(30) 



where P^t is the parallel-propagator from 71 (ri) to 72(1") along 
p(s, t). In the previous section, I did not write P'^b explicitly, but 
here it is best to do so, to guard against the temptation to use a 
background Minkowski structure to subtract (dri/dr)7i from 72 
(and thus neglect a potential first-order part). In general, I will not 
write the parallel-propagator factors in terms which are already 
first-order (since the omitted corrections would be of higher order), 
but I will keep them in zeroth-order terms. 
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For later reference, note that 



72 • V^" = 



-72 



dr2 



Pa • b 



dn . a . b 

——72 • VF i,7i 



= 72 



dri 
d7 



Pa ■• b 
b7i 



w''/''i?P,6"7i'ds' ) . (31) 



(In keeping with the remarks above, since the curvature is first- 
order, I have omitted the parallel-propagator terms which should 
properly appear in the integrand.) From this, we have 

r • v7^£i dri ,. r,a -b f dTl\ . -b 



dr' 
dn 
" d7 



d Tl f . r,a -b dn 



RpqabW^ I'' ilil ds 



(32) 



In these equations, and in others that follow, the field w'^ appears 
within terms which are already first-order. In such terms, we need 
only the zeroth-order expression for w"; this is 



a, S2 - S dn .a , S - Si 

V (s) = — 7i + —72 



(33) 



S2 — Si dr ■ S2 — Si 

which linearly interpolates between the values of ly" at the ends. 
(Because this will be used only when multiplied by first-order fac- 
tors, we do not write the parallel propagators which transport the 
vectors at the ends of the null geodesic to p{s, r).) 

The remaining quantities we shall need are the solutions [/"i,, 
V^b to the Jacobi equation, the vector /3" which is one of the ini- 
tial data (the vector a" — {dn/dT)^f is already known), and the 
luminosity distance r. 

To first order in the metric perturbation, we have 

U\ = {5\ + u\)P\ (34) 
V\ = ({s- s-,)5\ + v\)P\ (35) 

where now P'^b = P'^tis) is the parallel-propagator along p{s, r) 
from 7i('ri(r)) = p(si, r) to p{s, r) and 



u b = 



V b 



(s - s)H{s ~ syi^Rpiq^-ds (36) 
' (s - s)(s - si)H{s - s)Fl''Rpbq'"ds , (37) 



with H the Heaviside step-function, and the curvature is evaluated 
at p(s, r) in the integrands. 

We find, in the linear approximation, that 



/3" 



-{P ) bO u i,7i 



S2 - Si 



S2 - Si 



(S2 - Si)' 



a cb 
;V bO , 



(38) 



where on the right the first term contains the zeroth-order contribu- 
tion, and u";,, y'^b are evaluated at s = S2. 
Finally, from eqs. l[37j and p.0) , we find 



= 71 • ^ (^{s2 — Si) — v'^bmarn^ 
= 71 ■ ^ ^ (S2 — Si) 



(39) 



-AtyG 



where Tab is the stress-energy and we have used Einstein's equa- 
tion. 

3.1 Basic results and discussion 

With the results of the beginning of this section and a bit of work 
using eqs. jl l[ l, ^12\ , l |14[ l, we have 

dVi 
d72" 



(/•7i)"' 



(S2 — Sl) ^S-S+ / Fl''RpaqbV''v'' ds 



+ {l ■ jl) ■ ^2 ~ {I ■ jl) • 72)^^ ■ 72 ■ 



(40) 



A similar computation, using eq. ( |16^ , gives us the rate of 
change of the wave- vector fc" = ujiF. We have 



72 ■ Vfc" 



LJl (5" ^ CJl 

^ ■ 71 S2 — Sl r 71 



Ffv'' —Rvba"^ ds 



S2 - Sl 



/ • 72/ ■ 71 + 



I ■ 71 
-71c , 

S2 - Sl 

And for the change in field, we find 



^iciP-^Vbd" 



,p,q b ^2 — S 



S2 - Sl 

Rpbq'^ ds 



(41) 



• V7 r S2 - Sl / dn ,. 71-5 

72 • VFab = -3—71 • ^ H 

r V dr S2 — Sl 



7P C 7(7 . d ^2 S 

rv rji — 

S2 - Sl 



Rpcqd ds ] Fa 



AtvG 



+ 



dn 
d^ 



(S2 - S)(S -Sl)v VFl^Tpq ds ) Fab 

71 • V(rF,6)| 



-2 



Fv''Rpq[a'' ds ) Fh]c 



(42) 



Equations l |40^ , l[4TJ and lj42j describe basic observables, in 
the limit of linearised gravity. There are several points to make 
about these results: 

(a) The possibility of ascribing an observation of one the left- 
hand quantities to gravitational radiation relies on having some sort 
of extra information allowing one to discriminate between the var- 
ious terms on the right. In most cases, we must assume that time- 
dependence of the gravitational waves is enough different from 
those of the other physical processes that this can be used. 

(b) The equations contain terms proportional to 5" or 5 • ^ di- 
vided by i ■ 7i(s2 — Sl). These terms typically are suppressed when 
the observer is at very great distances from the emitter. (In fact, in 
view of the formulas l |31| ), l |32^ , \i9\ , for monochromatic waves of 
angular frequency uj, the first-order parts of these terms are typi- 
cally suppressed by factors of {ujt)~^.) In such cases, the zeroth- 
order contribution to d'^n/dr^ may become effectively negligible, 
leaving a geometrically pure first-order curvature -integral term. 

(c) In vacuum, the second term on the right-hand side of eq. 
^42\ will vanish. Note that, apart from this term, all of the curvature- 
integral terms are sums of 



r(") _ 

^ab ^ 



S2 - Sl 



/ / Rpaqb ds 



(43) 



for n — 0,1, 2, where we have used the formula \33\ for v°'. (Note 
that the form (Tsl of Fab imphes this for ((42l.) 
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(d) Many traditional approaches to these problems aim to 
work out what we might call the long-range scattering, correspond- 
ing to the receiver and emitter receding to great distances along the 
same null geodesicp^The formulas here show, however, that only 
in restricted circumstances will this limit exist. 

For the long-range scattering, we want to examine what hap- 
pens as si — >■ — CX3, S2 — >■ +00. Of course, we must assume that 
the contributions from the accelerations 71 , 72 are negligible (or 
at least are stable under the limit), and that the 5" -dependent terms 
drop out. But even then the remaining curvature integrals will not 
in general stabilise. This is because they depend on the quantity 
v" (s) (eq. j33| ) which interpolates from (dri /dr)7i at 71 to 72 
at 72. This quantity has no well-defined limit point- wise in s as 
si — > —00, S2 — >■ +00 independently. In other words, the contri- 
butions of the curvature integrals are in principle sensitive to the 
choices of of si and S2 in the asymptotic regime. 

On the other hand, in many cases of interest 7" and -yf will 
differ by only sub-relativistic effects, and then v"" (s) will be nearly 
constant along the null geodesic. Then the curvature integrals con- 
tributing to the time-dilationl|40[( and field change ^42\ will (as- 
suming the curvature falls off suitably) stabilise as si — > —00, 

S2 — >■ +00. 

Below, we shall mostly be interested in the case where the 
receiver is removed to arbitrarily great distances, but the emitter is 
held fixed. In this case we will have — >■ (dri/dr)7i, and only 
the integrals l[43j for n = will contribute. 

(e) Besides the implicit dependence of v'^{s) on si, S2, the 
curvature integrals in the expression (j4TJ for the change in wave- 
vector involve explicit factors (s — si)/(s2 — si), (s2 — s)/(s2 — 
si). These factors are of geometric origin, and express the fact that 
change in angle perceived by a distant observer will be of order half 
the full scattering angle multiplied by the ratio (distance of source 
to scatterer)/ (distance of source to receiver). In practice this means 
that if the source of the light is much closer to the source of the 
gravitational waves than it is to the Earth, the angular change due 
to the scattering is correspondingly reduced. 

Thus attempts to measure angular deflections due to gravita- 
tional waves are at a geometric disadvantage relative to measure- 
ments of changes in time-dilation or field. (This point can also be 
deduced from the formulas in |Kopeikin et al.|1999| ) 

(f) Below, we shall be interested in the case where the receiver 
is very distant but the emitter is not, so S2 — > +00 but si is finite. 
In this case we will have — > (dri/dT)7i, and the change in 
wave-vector will be suppressed. 



4 QUADRUPOLE SOURCES 

With the formulas derived above, the analysis of the effects at the 
linearised level in any given space-time reduces to the computation 
of certain moments of the curvature over the relevant segment of 
the light-ray's trajectory. The curvature can itself be expressed as 
a retarded field due to sources, plus a possible pure radiation term. 
The results of this can be quite complicated, even in simple cases, 
because of the time-dependence of the curvature and the different 
components which enter, and the fact that we wish to take the point 
of emission of the light to be finite. However, for the remainder of 



this paper, the aim will not be detailed modelling but simply rough 
estimates of the scales of the effects. 

I shall here work out the leading contributions in a simple but 
important case: a pure quadrupole field, from an isolated source, 
with the null geodesies in the radiation zone in the sense that 
ojb 2> 1, where uj is the angular frequency of any component of 
the gravitational wave and b is the null geodesic's impact parame- 
ter relative to the quadrupole source. 

The point of reception will be taken to be very far away from 
the source; this corresponds to the limit S2 — >■ +00 discussed ear- 
lier, but si will be held finite (recall si, S2 are the affine parameters 
specifying the null geodesic segment from emission to reception). 
In this case, the only curvature integral l |43| > we have to compute 
is /jI]'' (because within the integrands we have — >■ (dri/dr)7f , 

((S2 - S)/{S2 - Si)) ^ 1, ((s - Si)/(S2 - Si)) 0). 

Of course, real sources have monopole and perhaps dipole as 
well as quadrupole components; however, these contribute only sta- 
tionary terms to the field, and in any event in the linear approxima- 
tion those can simply be added to the quadrupole effects. 

4.1 Quadrupole fields 

By a quadrupole field I mean a linearised gravitational field in an 
appropriate j — 2 representation of the rotation group; both 'elec- 
tric' and 'magnetic' quadrupoles (often called mass quadrupole and 
current quadrupole terms) are allowed. The treatment here is cho- 
sen to fit with rest of this paper's formalism; other forms are given 
in |Regge & Wheeler|(T957) ; |Pirani|p965) ; |Thorne|(T980l l. 

Let us first consider the 'electric' part, which we idealise as a 
pure quadrupole at the spatial origin. (Since we are only interested 
in the field outside the source, this is adequate.) Let the quadrupole 
moment be Qa'(,(f), with arbitrary time-dependence. (Here t is the 
the coordinate time at the spatial origin, and is symmetric, 
trace-free and orthogonal to f".) It makes a contribution 

t:1 = {l/2){t-VS''a^taV'')it-V5\~UV'')Q;\(t)S^''^ (44) 

to the stress-energy, where 5'^' is the spatial delta-function. One 
easily verifies that / {tPt''Tpl){xa - tta){xb ~ ttb) d^x = Qf,, - 
the mass quadrupole is indeed Q^'i,. 

The full curvature tensor can easily be worked out by standard 
means [^It is convenient to introduce a null tetrad i", m", m", n", 
with = i" + f", for a unit space-like radially outward vector, 
=r-f'',andm" = 2"^''^(as -i esc 619^). Then the radiative 
(order r~^) term is 

Kbcd = -4G ^ Z[„m6]/[cmd] + conjugate +• ■ ■ , (45) 

where the superscript (4) indicates the order of differentiation with 
respect to u. (Note that the polarisation factors l[amb] are the same 
as in the electromagnetic case.) One does not actually need the de- 
tailed form of the m" vectors in computations; the combination 
mpTnt entering here may be written as 

nipmt = (l/2)(-gpi, + tptt — fpfb + iepbqst''f") . (46) 

Weyl curvatures of 'magnetic' type can be obtained, in lin- 
earised theory in the vacuum, by dualizing the electric ones. Thus 



" It might be tempting to call this the total scattering, but that would be 
misleading in this context, because we still track here only the differential 
effects, as the light-rays vary, of the scattering. 



One point to be careful of in these calculations is that /(u)Va5^^' 7^ 
/(t) Va^f''' in general, as becomes clear by multiplying by a test function 
and integrating by parts. 
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a 'magnetic' quadrupole field will, in the vacuum region, be given 
by 



mag (4) n a 



l\ami,]lyc.md\ + conjugate + ■ 



(47) 

where Q^g^ is referred to as the magnetic part of the quadrupole 
momentr^ The form of the source for this term is different from 
l|44j; one can check that the corresponding contribution to the 
stress-energy is 

(48) 

It turns out that Q^^^ is essentially a first spatial moment of 
the angular momentum density. To see this, note that = 
—epqr'^{if){x''){t'^TJ) can be interpreted as the angular momen- 
tum density with respect to the spatial origin (the minus sign giving 
the usual convention for the angular momentum as a spatial vector). 
Then a short calculation shows 



(49) 



One may take a complex quadrupole moment Qah = Qai + 
'QoiT^' the curvature in the radiation zone is 



Rated — — 4G 



/[ambj/jcmd] + conjugate H . (50) 



4.2 The light-rays 



We wish to study the differential scattering of light-rays which pass 
through the gravitational wave source's radiation zone and are re- 
ceived at some great distance. In this subsection, we work out the 
appropriate parametrization of those rays. Since we are to com- 
pute an integral of the curvature, which is a first-order quantity, it 
is enough to know the geometry of the ray and of the receiver to 
zeroth order 

Let the Bondi coordinates, centred at the world-line of the 
gravitational-wave source, be [u,r,9,<j)). Actually, we will not 
need to write [6, (j)) explicitly; we may represent them by their cor- 
responding null vector Z° = l°'(d, (f>), normalised hy I ■ t = 1, with 
the unit future-directed time-like vector characterising the Bondi 
frame. A point in Minkowski space is thus specified as ut"' + rZ". 
We will suppose the light-ray is received at an event 

72 (r) = U2t°- + r2l2 (51) 

where r2 is very large. 

In general, the equation of a light-ray may be expressed con- 
veniently as 

p''(s) = {u^h)e + bl'^ + sl1, (52) 

where u is the retarded time the light-ray tends to as s — >■ +oo, the 
null future-directed vectors /S- are normalised by Zo ■ = ■ ^ = 
Zi ■ f = 1 (so the spatial parts of ZJ}, are orthogonal); then b is the 
ray's impact parameter, its point of closest approach to the spatial 
origin occurs at s = 0, and Iq, code the direction of closest 
approach and the direction of the ray. We have r = + and 
r = (6/r)/S + {s/r)ll + (1 - 6/r - s/r)e. 



^■^ The sign here is fixed by the convention that it is the real and imaginary 
parts of the Bondi shear which determine the electric and magnetic parts of 
the curvature. 



light-emitter 




to observer 



gravitational-wave source 



Figure 1. The geometiy of the gravitational-wave source and the light- 
emitter The impact parameter is b. The coordinate along the light-ray is s, 
with s = corresponding to the ray's closest approach to the gravitational- 
wave source, and s\ the value at which the light-emitter sits; the observer is 
at a very large value S2 ■ The distance between the gravitational-wave source 
and a point on the light-ray is r, with r\ the distance of the light-emitter 
Then angle 9 (not used until eq. j76}) is given by si = ri sin 9; note that 
s\ and 9 may have either sign. 



Requiring the ray to be received at p"(s2) = 72 (''")' we find 

u + S2 = U2 + r2 (53) 
h^ + sl = rl (54) 
[b" + sl)-^/\bro+S2ll) = /a" modi", (55) 

where 'mod f' means up to terms proportional to t". Solving these 
equations perturbatively in 6/r2, we have 

S2 = r2 + 0{b^/r2) (56) 
u = U2 + 0(S?/r2) (57) 
n = /^ + 0(6/r). (58) 



To this order, the vector Iq is (apart from the normalisations speci- 
fied above) unrestricted; the quantity {b/r2)lo specifies the appar- 
ent direction of the light-ray, relative to the direction of the source, 
at the receiver. 

The light-ray's trajectory is thus 



— ('"2 — 6)i" + Uq + SI2 , 



(59) 



meeting the receiver at s = S2 = '■2- The point of emission will be 
at a parameter value si, the point of closest approach to the spatial 
origin, as noted above, would be s = 0. Expressing p''(s) in Bondi 
coordinates, we have 

p"(s) ^Mf^+r/", (60) 
where 

u = U2 + S- v/b^Ts^ (61) 
r = \/62 -f s2 (62) 
= {b/r)ll + {s/r)l1 + {l-b/r- s/r)e . (63) 

In the computations to follow, it will be convenient to put 
s = 6 sinh ^ . (64) 
Then r — 6 cosh ^ and bs — y/h^ + — —be~^ and we define 

S = -e-« = (s/b) - yiTTW ; (65) 
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note that S is an increasing function of s, and that 

ds — b cosh ^ 

_ _ dS^ 
4.3 The curvature integral 

As noted above, the only curvature integral we require is 



(66) 



r(0) _ 



(67) 



where = l\- The integrand, in the radiation zone, is 



l°'l'^Rab 



Ar^^pq laf ' -^v,c{ - -^1 
— 4G — —L L^am^m I l[^m^]m 

+ conjugate + ■ • ■ , 



(68) 



where the accents indicate the vectors evaluated at s along the null 
geodesic. 

Because we wish to evaluate J^|]' in the case cjb ^ 1, where 
to is the angular frequency of any contributing component to the 
gravitational radiation field, we shall for the moment just work with 
one Fourier component, putting e^'^^Kpq (where Kpq is constant, 
symmetric, trace-free, and orthogonal to f") in place of Qpq', after 
using the condition cjfo ^ 1 we will restore Qpq. Then making use 
of eqs. l|61[) and \65\, we find 



r(0) 



4Goj^ / e 

Jsi 
+ conjugate 



-iu:(u2-\-bS) 



K pql I 



S 
(69) 



In this form, the integral is proportional to e^"^ , and we may 
regard it as effecting a Fourier transform. We are interested in the 
behaviour of this for large ub, which is to say the high-frequency 
regime. The function to be transformed is smooth except for being 
cut off at the end-points 5 = ^i, 5 = 0; it is the non-smooth 
behaviour at these end-points which will give the leading contribu- 
tion. 

In fact, it is the lower end-point S = Si which makes the 
dominant contribution, for a little algebra shows the integrand tends 
continuously to zero as S ^ 0- We have then the Fourier transform 
of a function with a jump discontinuity at S — Si; this is 



r(0) 



"UGcj'^Kp 



f 

Js, 



-iu}(u2'\-bS) ja 



,a, , —Pjc{ - -JL.q dS \ 

I liairiMm I L^m^^m — 

' ' ^ ' S ls=5i 



+ conjugate + • • • 

-AGQ'fqriya'm^Tfi^ flycThifm^ {bSy 
+ conjugate + • ■ ■ , 



\S=Si 



(70) 



where the tilde denotes asymptotic expansion for ub ^ 1 and we 
have restored Qpq . 

We will be most interested in the case of time-delays, for 
which the curvature-integral contribution is 



d^ 



{dri/drfl^:jtH\ 



(71) 



Then the combination of vectors entering into eq. l |70^ becomes 

4t''l'^l^ami,]fn' = —2l'^rhafrf\ (72) 

I s — SI 

^ ( ^IP ^ i^C^M mod^^ 

r \r r J r \ ij=si 

where terms proportional to have been dropped because they will 
not contribute when contracted with Qpq and 



O — fc aqsi'l C to 



(73) 



is a purely spatial vector equal to lo x li, where lo, li are the spatial 
parts of Writing now ri = r = ^/b^ + (^i)'^, we have 



d^i /dnV Gb^ 

dr'-^ \dr J 4(ri — Si)(r-i)2 

ri J \ri 
+ conjugate + 



X Qf^ ( 



,q _ Si q 
«1 «0 



ri 



(74) 



or 
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b q _ Si q 

n ri 
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Wpq 



b ,p si p 

'l '0 

n n 



ri 



H 



■ (75) 



(3) 

Here Qpq is evaluated at it = U2 + Si — ri. 

We may recast the foregoing in terms of the angle 9 of the 
light-emission from the point of closest approach relative to the 
source (see Fig.[TJ, so 



Si = ri sin 9 , b = ri cos 9 

(where b is the impact parameter); then we have 



(76) 



dVi 
d72" ' 

1 

2' 



dri 
d7 



G(l + sin6i) 



ri 



)cl(3) 
,pq 



{l^ COS 9 — Iq sin 9) {ll cos t 



Iq sin ^ 



cos" 



'Q 



mag (3) 

pq 



{ll cos 9-11 sin Q) {C cos I + • ■ ■ . (77) 

Formula \n\ is the main result, relating the clock times ri of 
the emitter (at radius ri) and r of the receiver, as influenced by 
gravitational mass ('electric') Q^';, and current ('magnetic') Q^jf^ 
quadrupole sources (evaluated at retarded time u = U2 + si—ri — 
1*2 — (1 — sin&)ri, where U2 is the observer's retarded time). Re- 
call that here li is the null geodesic's tangent and Iq is the null 
vector whose spatial part lo is a unit vector from the source to 
the geodesic's point of closest approach, and C" is a spatial vector 
Iq X ll normal to the plane containing the source and the light-ray. 

Perhaps the most striking feature of this result is the 'forward- 
backward' asymmetry represented by the overall factor (1 + sin 9), 
which enhances effects from light-emitters on the portion of the 
light-ray outgoing from the gravitational-wave source (0 < ^ < 
7r/2) relative to those from the incoming portion (— 7r/2 < ^ < 0). 
This is a relativistic effect arising from the use of light-signals to 
probe the space-time curvature. For gravitational waves of a given 
frequency with respect to (the gravitational source's frame). 
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the frequency with respect to an affine parameter along the hght- 
ray will be larger along the ingoing portion than along the outgo- 
ing one, and the effects due to those higher-frequency terms more 
nearly average out. 

For the scaling of eq. l |77| l with distance, for light-emitters 
very distant from the light-ray's point of closest approach to the 
gravitational-wave source, that is ri ^ 6, we have 1 + sin 6 = 
1 ± y/1 - {b/rxY and 



1 + sin t 
ri 



h^/{2rl) for6»i-7r/2 
2/ri for6ltvr/2. 



(78) 



Thus the rough magnitude of the effect will be 



'|d-ri/rfr|^ 




{GbVirfc^)) 


fore i -tt/2 






(G/(bc^)) 


for moderate ( 


^ 1 






fore* t vr/2, 



dr2 



(79) 

where the speed of light has been given explicitly. (For rough esti- 
mates, the precise choice of norm for the tensor is not very impor- 
tant; any Lp norm in terms of a standard Euclidean basis will do.) 
While the upper line corresponds to the scaling found by Damour 
and Esposito-Farese, we see that the fall-off for light-sources along 
the outgoing portion of the ray is much softer, having the ~ 
behaviour characteristic of radiative effects. 

While the appearance of this radiative scaling is certainly of 
interest, we shall see below that even in favourable circumstances 
the effects are small; we shall therefore concentrate, in the follow- 
ing sections, with the case of moderate 6, corresponding to light- 
emitters with ri ~ b, the middle line of eq. \79\ . 

Finally, two remarks about the angular dependence of the 
effects through the factor in curly braces in eq. |77]l. First, the 
variations of this term for moderate 6 means that different light- 
emitters in this regime will probe the different components of the 
quadrupole tensors. Second, one might have thought that, for light- 
sources further away, the trigonometric factors would tend to sup- 
press the dependences on Iq, the vector from the origin to the point 
of closest approach, and lead to a dependence of the effects primar- 
ily on li, the tangent to the light-ray. The opposite is true of the 
factor in curly braces in eq. (|77j, however. This is a direct conse- 
quence of the trans versality of the waves; the components of the 
curvature that enter are orthogonal to the position-vector relative to 
the origin. 



5 ESTIMATES 

Of the various possible modulations of the light by its passage 
through, and emission within, the gravitational radiation - changes 
in the received light's amplitude, polarisation, phase, location on 
the receiver's sky, and time-dilation - it seems that in most cases 
time-dilation will be the most promising for detection (but still, as 
we shall see, quite challenging). As pointed out earlier, changes 
in the location and phase will be suppressed by a geometric fac- 
tor. Changes in amplitude and polarisation would be so small they 
would probably be too hard to detect. 

In this section I shall estimate two time-dilation effects, red- 
shift and pulse time-of-arrival offsets, in some cases of interest. 
First some rough general formulas for these will be derived; we will 
see that they are expressed naturally in terms of a dimensionless 



intrinsic measure of gravitational-wave strength, the Bondi news 
(essentially GQf^ /c"' in our case). While this quantity is central 
in much of gravitational radiation theory, few numerical values for 
it have appeared in the astrophysical literature, so samples of these 
are given. 

I then discuss the pulsar timing resolutions appropriate to the 
detection of time-of-arrival wanderings (unfortunately, the extraor- 
dinary long-term stability of pulsars does not help directly with 
this), and finally estimate those wanderings. 

5.1 Formulas for the red-shift and time-offsets 

The entire light-signal from the emitter is subject to a time-dilation, 
which will itself be a function of time. If the emitted signal is a(ri) 
in the frame of the emitter, then, the received signal will be (ne- 
glecting other effects) a(ri(r)); that is, there will be a distortion 
due to the relative differences in the flows of time. One effect this 
would give rise to would be a time-dependent red-shift. Or if one 
knew the light-source was, in its own frame, emitting regular sig- 
nals (for instance, if it were a pulsar), the effect of the gravitational 
radiation would be to make the times of receipt of these wander 
slightly from complete regularity. 

Where linearised gravity is adequate, the equation governing 
this was \4Q\ , and it contained three sorts of contributions: terms 
due to possible covariant accelerations of the emitter or the re- 
ceiver; a sort of kinematic term due to the possible boost of the 
receiver relative to the emitter; and the curvature integral estimated 
in the previous section. It is the last which will be important, as will 
now be explained. 

Recall that our main hope for detecting gravitational waves by 
this effect comes not from relative magnitudes of these terms (the 
gravitational acceleration at the surface of the Earth is far larger 
than the expected gravitational-radiation effects), but from the dif- 
ferent time-dependences of the terms. We must assume that we can 
account for any source and receiver accelerations well enough to 
distinguish the effects from those of possible gravitational waves. 

The kinematic term (5-5/(/-7i(s2 — si)) requires a bit more 
discussion, though. Here 5° is the difference between (dri / dT)jf 
and 72, parallel-propagated along the null geodesic. While in the 
situations we shall consider the zeroth-order contribution to this 
will have a different time-dependence than the gravitational-wave 
effects, the first-order contribution will be affected by the gravita- 
tional radiation. However the corresponding effects are suppressed 
by a factor of (cjr/c)"^ (essentially because they involve velocities 
rather than accelerations), as inspection of eqs. l |32[ l, (j4TJ shows; 
see the discussion under point (b) in 3.1. 

We thus consider only the curvature-integral contributions to 
the time-dilation. Assuming that the light-source is not moving 
ultra-relativistically with respect to the Earth, we see from eq. l |79| l 
that the scale of the effects is, in the case of light-sources ~ b from 
the gravitational-wave source, where b is the light-ray's impact pa- 
rameter. 



d72" 



- Ilo'^'l 
b V c5 ir"" \ 



(80) 



where the first factor on the right has units inverse time and the 
second is dimensionless (recall is the third time-derivative of 
the quadrupole moment) 

The second factor on the right in eq. JSOl is a dimensionless 



The right-hand side of eq. jSOl also gives a rough estimate of the frac- 
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measure of the intrinsic strength of the gravitational radiation; it is 
essentially an average value of the magnitude of the Bondi news 
N, a key quantity in gravitational radiation theory which will be 
discussed a bit further below. Because this is such central concept, 
we now switch to expressing the quantities of interest in terms of 
this average news 

G 



\NU 



„5 \\^P1 



Then the rough magnitudes of the red-shifts will be 



c 

Lob 



\NU 



and the magnitude the time-of-arrival wandering will be 



I At 



c 



\Nl 



(81) 



(82) 



(83) 



where uj is the angular frequency of the wave and b is impact pa- 
rameter of the null geodesic from the gravitational wave source. 
(Again we assume the emission occurs near the point of closest ap- 
proach - more generally one would replace the factors of 1/6 with 
(1 + sin6)/r-i, in accord with \n\ and that the motion of the 
receiver relative to the emitter is sub-relativistic.) 

We recall that the analysis of the previous section assumed 
that the light emission occurred in the gravitational wave zone. 
This means that we must have ujb/c > 1. (Emissions from points 
closer to the gravitational source could be studied by the gen- 
eral formulas given earlier, but they would correspond to near- or 
intermediate-zone gravitational disturbances, which had not prop- 
agated far enough for their wave character to be fully developed.) 
Thus we have 



I 2 1 max 

and 

lArlma 



(84) 



(85) 



for rough estimates of the largest possible red-shift and time-offset 
due to light emissions near the geodesic's point of closest ap- 
proach, in the gravitational wave zone, for a given gravitational- 
wave source. 



5.2 Source types and intrinsic wave strengths 

In much of the literature, gravitational radiation is estimated by a 
combination h of the linearised metric components at the Earth. 
While useful for discussions involving terrestrial detectors, this 
measure is neither invariant nor intrinsic; it is not well-suited for 
the present considerations. 

The measure which is appropriate is the Bondi news N, which 
we have already mentioned. This is a function of retarded time and 
angle; it has a suitable invariance (Bondi-Metzner-Sachs covari- 
ance), is intrinsic and dimensionless (Pe nrose & Rindlerl|1986[ l. 
The gravitational luminosity is (c^ /4ttG) § \ N\ (over the sphere 
of directions). In the quadrupole approximation one has = 
— {G/c^)Q^^,jm''fn'' (where m" is a complex null vector coding 
the direction), and so the quantity jA'^jav given in eq. jsTJ is an av- 
erage of lA'^l over directions. There are, however, few numerical 
values of the news in the literature. 

Table 1 gives some rough estimates for the news, using 



tional rate of change of the electromagnetic field components, induced by 
the light's passage through, and emission in, the gravitational radiation. 



i A'^lav = (GL/c^)^''^ where L is the gravitational luminosity (and, 
for the last four lines, the quadrupole approximation). Note that the 
first three cases correspond to burst-type sources. 

It is clear from Table 1 that if we should be lucky enough to 
observe light-signals from the vicinity of an extreme gravitational- 
wave event, such as an asymmetric supernova or colliding black 
holes of comparable mass, a great deal of information could 
be gained. However, because the gravitational waves decay very 
quickly in such cases (typically on a time-scale of order the light- 
crossing time associated with the mass, that is GM/c^), one would 
have to already have the telescope trained on the object, and in 
general this would require extraordinary serendipity. An exception 
would be if we could detect an inspiralling system and thus be pre- 
pared to monitor light-sources in its vicinity when strong gravita- 
tional waves were produced. 

For the rest of this paper, we leave aside the possibility of 
gravitational-wave sources with news any significant fraction of 
unity. The best other candidates for detection appear to have bi- 
naries as sources of gravitational waves, corresponding to the last 
two lines in Table 1 . For these, we have 



liVL 



1.3 X 



fGM_ 



27r 



5/3 



for the case of equal masses in a circular orbit and 



\N\, 



27rGM 



2/3 



P I 



(86) 



(87) 



for highly unequal masses M 3> m in a circular orbit; in each case 
P is the period. 

5.3 Red-shifts 

For neutron-star glitches, the factor 8 x 10^^ in the third line of 
Table 1 corresponds to a spectral resolving power of 20kms^^. 
The red-shift will be smaller by a factor (c/ujh), and for a light- 
source ten light-seconds away from the glitching neutron star (say, 
a companion) this is likely ~ lO"'' (taking l/o; ~ 1 ms), giving 
\z\ ~ 2ms^^/c. Radial components of st ellar velocities are mea- 
sured, interferometrically, tol — 3ms~^ (Butler et al. 1996 Rup- 



Iprecht et al.|2004| l; however, such precisions typically require inte- 
gration times of a minute or more, whereas the glitches are thought 
to decay over milliseconds. 

For the other sources we considered, the red-shifts would be 
smaller, and hence we now move on to the possibilities of detecting 
the wanderings of times-of-arrival of signals from pulsars. 

5.4 Pulse time-of-arrival offsets 

We now come to the possibility of measuring the offsets caused 
in pulsar signals times-of-arrival caused by their passage through 
gravitational waves. We saw above (eq. \%5) ) that the maximum 
offset was roughly Lj^^|Af|av. The effect is therefore determined 
by a ratio of the intrinsic gravitational wave strength |Af|av to the 
gravitational wave period oj, and we must consider this for different 
classes of sources. 

First, however, let us consider what temporal resolution we 
might get from the pulsar signals. 

5.4.1 Millisecond pulsars as clocks 



As account of pulsar and timing is given by |Lorimer & Kramer| 
(j2005j. While millisecond pulsars can have extraordinary long- 
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Table 1. Rough estimates of the Bondi news |A'^|av for some classes of gravitational-wave sources, inferred from estimates 
of the gravitational luminosity (and, for the last four classes, the quadrupole approximation). The first three classes are burst 



Source type 



Estimate 



Parameters and expected ranges 



highly asymmetric supernovae 
mergers of comparable-mass black holes 

neutron-star glitchei 



long-duration waves from neutron stari 

binary, equal masses At, circular orbit 
binary, masses M 2> m, circular orbit 







possibly a significant fraction of unity 
possibly a significant fraction of unity 
p \ 1/2 



X 10" 



^ 10-'' Mqc- 



/ 1ms \ 



1/2 



1 X 10" 



. looAi-Q ; 



Id 



5/3 



Egvi = energy in the wave 

Tdccay = decay time 
IO-^Mq < Egw < IO-^Mq 

e = efficiency 

e < 2 X 10^^ for periods as long as years 
P = period 
P = period 



Adapted fromlAndersson et al.'l 20 1 1 1. 

Adapted from [Kokkotas et al.|^2001^ , |Prix 



2009 



term stability, that stability is essentially a secular effect, based on 
measuring the times and numbers of pulses between initial and final 
pulses years apart, and having only to reckon with the uncertainties 
in locating the times of arrival of those initial and final signals. Ex- 
cept in searching for very low frequency gravitational waves (or 
'memory' effects, which are essentially zero-frequency phenom- 
ena), that stability is not of direct help. We must consider timing 
over shorter scales. 

There are two key issues: that pulsars are weak radio sources, 
and that their signals over a period appear to be subject to in- 
trinsic variations, jitter and drifting sub-pulses. Because of these, 
one cannot extract very accurate times of arrival from individual 
pulses; instead one must fold them, that is, integrate (typically 
over a few minutes) with a mean periodicity removed in order to 
get a good pulse profile. One may then have a resolution around 

If we search for effects corresponding to shifts in the times of 
arrival of order At Arpuisar, we need of order (Arpuisai/ Ar)^ 
samplings, each of some minutes' length. There are about 5 x 10^ 
minutes in a year, and so it would be very hard to track changes 
smaller than ~ lO"^*^ s by this crude statistical means with current 
technology. 

However, were we to learn more about the jitter and sub-pulse 
drift, we might be able to do better. Most of the effects consid- 
ered here will have fairly well-defined periods, and the question is 
whether those Fourier components could be distinguished from the 
jitter and drift. These displacements of the times-of-arrival would 
show up in the Fourier domain as splittings of the angular frequen- 
cies of the components of the signal by ±ix), or as peaks in the 
Fourier-transformed residuals, at angular frequencies ±cj. 



5.4.2 Long-duration waves from neutron stars 

The case where the pulsar is one member of a binary, the other 
member being a neutron star gravitationally radiating due to an 
asymmetry or instability (fourth line of Table 1), turns out not to 
be very promising. (Note that here the gravitational radiation in 
question is not due to the binary system.) Taking the case of a non- 




10 

M (solar masses) 



100 



Figure 2. Contours of approximate maximum time-offsets Armax from 
two equal masses, each M, in a circular orbit of period P. 



axisymmetric neutron star considered by |Prix| ( [2009^ , we have 

At ~ liV|av(c/a;'6) 



(6 X lO^S) |iV| 



10 ms 



lOlt-s 



(88) 



where P is the period of rotation. Using Table 1, we see that in op- 
timistic circumstances for persistent waves, we might have effects 
as large as ~ 4 x 10"^® s. This is very small. 



5.4.3 Pulsar detection of binary emission 

It is more promising to consider as a source of gravitational waves 
a tight binary, and assume that there is a nearby pulsar. This pulsar 
could be a (relatively distant) tertiary, or not be directly gravitation- 
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Figure 3. Contours of approximate maximum time-offsets Armax from 
two unequal masses, with primary mass _A/ 2> m and secondary inass 
m = Mq, in a circular orbit of period P. 

ally bound but still close by (as might happen if, for example, both 
the binary and the pulsar were in a globular cluster). Recalling that 
cjb/c > 1 for the wave character of the binary-induced changes in 
curvature to be developed at the pulsar, and using eqs. ([86j, l |85[ l, 
we see that for a binary of equal masses, each A/, in a circular orbit, 
the maximum offset in pulse time-of-arrival will be 

A.„„L_-.(3xio-".)(^)°"(i^)". m 

Figure |2] shows contours for Armax for l |89[ l. For example, if 
we had the component masses M = SM© and period P = .Id, 
we would have Armax ~ 2 x 10~^°s. This (rough) estimate 
of the maximum time-displacement corresponds to the pulsar be- 
ing at the inner part of the gravitational radiation zone, that is 
~ cP/ (27r) ~ 30 1-min away from the binary. If the pulsar were 
as far away as 10~^ pc (representative of interstellar distances in a 
globular cluster), one would have Ar ~ 3 x 10"^"^ s. 

In the case M 3> m of one mass dominating the other we 
have (from eqs. l[87j, l[85j) 

(90) 

Contours of Armax for this case, for the secondary of mass m — 
Mq, are given in Fig. [3] For instance, if one could find a black 
hole of mass M = 10* M© (in, say, a globular cluster), with a star 
with mass m = Mq orbiting in a one-day period, one would have 
Armax ~ 3 X 10~^ s; again, this limit is set when the pulsar is 
at the inner part of the gravitational wave zone, in this case a few 
light-days away from the radiating binary. For a solar-mass star in a 
ten-year circular orbit about the M — 4x 10® A/q presumed black 
hole at the galactic centre, we should have Armax ~ 7 x lO"^'' s, 
and this would apply to pulsars around cP/{2n) ~ 1.6 1-yr from 
that system. 

It should be emphasized that these estimates are very rough. 



6 SUMMARY AND DISCUSSION 

The first main goal of this paper was to present an invariant frame- 
work for treating the differential scattering of light-rays in exact 
general relativity. This allows one to keep track of the different 
physical contributions to the various scattering effects, and to fo- 
cus on quantities of direct physical interest. 

The result of greatest near-term possible observational conse- 
quence is that light emitted from the vicinities of gravitational-wave 
sources may be scattered by much larger amounts than those dis- 
cussed by Damour and Esposito-Farese, and by Kopeikin, Schafer, 
Gwinn and Eubanks. The results here are however compatible with 
those authors'; one can understand them as due to 'edge effects', 
whose possibility was explicitly noted by Damour and Esposito- 
Farsee, and which could be treated within the general formalism of 
Kopeikin et al. An equivalent statement is that for light emissions 
at finite distances from the gravitational-wave sources, the scatter- 
ing depends on the gravitational radiation field, in contrast to the 
cancellations found by Damour, Kopeikin et al. which occur in the 
limit of infinite distances. The effects do come out to be, roughly, 
of the scale of those predicted in some of the still earlier papers 
l |Sazhin|I978[|Fakir|1994| l, and in a very rough sense one may say 
that this is because those papers did have the scattering respond 
to the gravitational radiation fields; however (again in accord with 
Damour, Kopeikin et al.) the details of the arguments of those pa- 
pers are not really compatible with those here. 

The main physical issue which was not addressed systemati- 
cally in earlier work was the effect of the gravitational waves on the 
motion of the light-sources. (See however [Kopeikin et al.|20rT| for 
an exception.) This showed up, mathematically, in leaving the re- 
sults in terms of the coordinate times rather than clock times of the 
light-emitter and receiver. Such formulas are of direct physical sig- 
nificance only if the coordinates are adapted to the motions of both 
the light-emitter and receiver. This issue is particularly problematic 
when the light-emitter is close enough to the gravitational-wave 
source that that gravitational effect on the light-emitter's motion 
must be accounted for. The present approach works with the clock 
times and can account for any motions of the light-emitter and re- 
ceiver. 

Some candidates for the possible observation of these scatter- 
ing effects were considered, in a quadrupole approximation for the 
gravitational waves. It was found that the best effects to search for 
were based on the relative distortion of the emitter's proper time n 
as measured by light-signals incoming to the receiver at its proper 
time r. 

Because these observations would require training a telescope 
on the emitter, one is driven to look for specific likely sources of 
gravitational radiation, and these are in general weak ones. The 
most promising sort identified here would be a tight binary, the 
light-emitter being a pulsar orbiting the binary as a distant ter- 
tiary, or not gravitationally bound but close by. In an example with 
favourable but not extreme parameters (a binary of two 5Mq stars 
with a . 1 d period, and a pulsar about 30 1-min away), we found 
offsets in the times of arrival of the pulsar pulses of the order of 
^ 10~^" s, with the offsets varying with a period around an hour. 
If we a were able to find a solar-mass star orbiting a 10* Af© black 
hole with a one-day period, for example in a globular cluster, pul- 
sars a few light-days away might have their pulse times-of-arrival 
offset by ~ 10"'^ s. 



© 2012 RAS, MNRAS OOO.pTfIS] 



Light, gravity and pulse-time offsets 15 



ACKNOWLEDGMENTS 

I thank Sergei Kopeikin and Bahram Mashhoon for drawing my at- 
tention to earlier work on this and for a helpful discussion. I thank 
Duncan Lorimer for a helpful electronic exchange about pulsar 
timing measurements, and an anonymous referee for constructive 
suggestions about the exposition. The graphs were prepared with 
Veusz. This work was supported in part by the University of Mis- 
souri Research Board. 



Sazhin M. V., 1978, SvA, 22, 36 

Schutz B. R, 2010, in S. A. Klioner, P. K. Seidelmann, & 
M. H. Soffel ed., lAU Symposium Vol. 261 of lAU Symposium, 
Astrometric and timing effects of gravitational waves, pp 234- 
239 

Shapiro 1. 1., 1964, Physical Review Letters, 13, 789 
Taylor J. H., Weisberg J. M., 1989, ApJ, 345, 434 
Thorne K. S., 1980, Reviews of Modem Physics, 52, 299 



REFERENCES 

Andersson N., Ferrari V., Jones D. I., Kokkotas K. D., Krishnan 

B., Read J. S., RezzoUa L., Zink B., 201 1, General Relativity and 

Gravitation, 43, 409 
Book L. G., Flanagan E. E., 2011, Phys. Rev. D, 83, 024024 
Bufler R. P, Marcy G. W., Williams E., McCarthy C, Dosanjh P, 

Vogt S. S., 1996, PASP 108, 500 
Damour T, Esposito-Farese G., 1998, Phys. Rev. D, 58, 044003 
Durrer R., 1994, Phys. Rev. Lett., 72, 3301 
Fakir R., 1994, Phys. Rev. D, 50, 3795 
Faraoni V., 2008, New Astronomy, 13, 178 
Hellings R. W., Downs G. S., 1983, ApJ, 265, L39 
Kaiser N., Jaffe A. H., 1997, ApJ, 484, 545 
Kokkotas K. D., Apostolatos T. A., Andersson N., 2001, MNRAS, 

320, 307 

Kopeikin S., Efroimsky M., Kaplan G., 2011, Relativistic Celes- 
tial Mechanics of the Solar System. Wiley. 

Kopeikin S., Mashhoon B., 2002, Phys. Rev. D, 65, 064025 

Kopeikin S. M., Schafer G., Gwinn C. R., Eubanks T. M., 1999, 
Phys. Rev. D, 59, 084023 

Labeyrie A., 1993, Astronomy and Astrophysics, 268, 823 

Lesovik G. B., Lebedev A. V., Mounutcharyan V., Martin T, 
2005, Phys. Rev. D, 71, 122001 

Lorimer D. R., Kramer M., 2005, Handbook of Pulsar Astronomy. 
Cambridge University Press 

Misner C. W., Thorne K. S., Wheeler J. A., 1973, Gravitation. 
W.H. Freeman and Co. 

Penrose R., Rindler W., 1984, Spinors and space-time, vol. 1: 
Two-spinor calculus and relativistic fields. Cambridge Univer- 
sity Press 

Penrose R., Rindler W., 1986, Spinors and space-time, vol. 2: 
spinor and twistor methods in space-time geometry. Cambridge 
University Press 

Pirani F. A. E., 1965, Introduction to Gravitational Radiation The- 
ory (Notes by J. J. Marek and the Lecturer). Prentice-Hall, p. 249 

Prasanna A. R., Mohanty S., 2002, FPL (Europhysics Letters), 60, 
651 

Prix R., 2004, in W. Beecker ed.. Astrophysics and Space Science 
Library, Vol. 357, Gravitational Waves from Spinning Neutron 
Stars, p. 65 1 

Regge T, Wheeler J. A., 1957, Phys. Rev., 108, 1063 
Rupprecht G., Pepe F, Mayor M., Queloz D., Bouchy F., Avila 
G., Benz W., Bertaux J.-L., Bonfils X., Dall T, Delabre B., 
Dekker H., Eckert W., Fleury M., Gilliotte A., Gojak D., Guz- 
man J. C, Kohler D., Lizon J.-L., Lo Curto G., Longinotti A., 
Lovis C, Megevand D., Pasquini L., Reyes J., Sivan J.-P, Sos- 
nowska D., Soto R., Udry S., Van Kesteren A., Weber L., Weilen- 
mann U., 2004, in A. F. M. Moorwood & M. lye ed., Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Se- 
ries Vol. 5492, The exoplanet hunter HARPS; performance and 
first results, pp 148-159 



© 2012 RAS, MNRAS 0Q0,[T|fl5l 



